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The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes 
I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations 
between the temperature and polarization maps significant at more than 10 standard devi- 
ations. The correlations are inconsistent with instrument noise and are significantly larger 
than the upper limits established for potential systematic errors. Correlations on small an- 
gualr scales are consistent with the the signal expected from adiabatic initial conditions. We 
detect excess power on large angular scales consistent with an early epoch of reionization. A 
model-independent fit to reionization optical depth yields results consistent with the best-fit 
ACDM model, with best fit value r = 0.17 ± 0.04 at 68% confidence, including systematic 
and foreground uncertainties. 

1 Introduction 

The Wilkinson Microwave Anisotropy Probe has mapped the full sky in the Stokes J, Q, 
and U para meters on angular scale s 9 > 0?2 in 5 frequency bands centered at 23, 33, 41, 61, 



and 94 GHz ([Bennett et al.Ll200 3a). WMAP was not designed solely as a polarimeter, in the 



sense that none of its detectors are sensitive only to polarization. Incident radiation in each 
differencing assembl y (DA) is split by an orthomode transdu cer (OMT) into two orthogonal 
linear polarizations ( Page et al. . 2003bt Jarosik et al. . 20031 ). Each OMT is oriented so that 
the electric field directions accepted in the output recta ngular wavegui d es lie a t ±45° with 
respect to the yz symmetry plane of the satellite (see iBennett et al.l (J2003aJ) Fig. 2 for 
the definition of the satellite coordinate system). The two orthogonal polarizations from 
the OMT are measured by two independent radiometers. Each radiometer differences the 
signal in the accepted polarization between two positions on the sky (the A and B beams), 
separated by ~ 140°. 

The signal from the sky in each direction h can be decomposed into the Stokes pa- 
rameters 

T(n) = 1(h) + Q(h) cos 2 7 + U(h) sin 2 7 , (1) 

where we define the angle 7 from a meridian through the Galactic poles to the projection 
on the sky of the E-plane of each output port of the OMT. Denoting the two radiometers 
by subscripts 1 and 2, the instantaneous outputs are 



ATi = I(h A ) + Q(h A ) cos 27^4 + U (h A ) sin 27^ 
- I{h B ) - Q{h B ) cos 27 B - U(n B ) sin 2^3 



(2) 



1 WMAP is the result of a partnership between Princeton University and NASA's Goddard Space Flight 
Center. Scientific guidance is provided by the WMAP Science Team. 
2 Code 685, Goddard Space Flight Center, Greenbelt, MD 20771 
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and 



The sum 



AT 2 = I{n A ) - Q{fi A ) cos 2^ A - U (h A ) sin 2^ A 
- I(n B ) + Q(n B ) cos 27 B + U(n B ) sin27 B . 

ATj = i(AT\ + AT 2 ) = I(n A ) - I(n B ) 



(3) 



is thus proportional to the unpolarized intensity, while the difference 



AT F 



-AT 2 ) = Q{h A ) cos 27a + U(n A ) sin 27^- cos 2^ B -U(n B ) sin 27^. (4) 



is proportional only to the polarization. We produce full-sky maps of the Stokes J, Q, and 
U parameters from the sum and difference time-ordered data using an iterative mapping 
algorithm. Since the polari zation is faint, the Q a nd U maps are dominated by instrument 



noise and converge rapidly (jHinshaw et all l2003a[ ) 



The Stokes Q and U components depend on a specific choice of coordinate system. 
For each pair of pixels, we define coordinate-independent quantities 



Q' = Qcos(2</>) + C/sin(20) 
U' = f/cos(20) -Qsin(20), 



(5) 



where the angle <fi rotates the coordinate system about the outward-directed normal vec- 
tor to put the merid i an al o ng the great circle c o nnect ing the two positions on the sky 



( Kamionkowski et all Il997t IZaldarriaga &: Seliaki [1997). All of our analyses use these 



coordinate-independent linear combinations of the Q and U sky maps. 



2 CORRELATION FUNCTION 

The simplest measure of temperature-polarization cross-correlation is the two-point angular 
correlation function 

C IQ {6) = 3 3 \ (6) 

where i and j are pixel indices and w are the weights. To avoid possible effects of 1/ f noise, 
we force the temperature map to come from a different frequency band than the polarization 
maps, and thus use the temperature map at 61 GHz (V band) for all correlations except the 
V-band polarization maps, which we correlate against the 41 GHz (Q band) temperature 
map. Since WMAP has a high signal-to-noise ratio measurement of the CMB temperature 
anisotropy, we use unit weight {wi = 1) for the temperature maps and noise weight (wj = 
Nj/ctq) for the polarization maps, where Nj is the effective number of observations in each 
pi xel 7 and an is the stan dard deviation of the white noise in the time-ordered data (Table 1 
of lBennett etaD (|2003bh . We compare the correlation functions to Monte Carlo simulations 



of a null model, which simulates the temperature anisotropy using the best-fit ACDM model 
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( Soergel et all l200aj) Dut forces the polarization signal to zero. Each realization generates 



a CMB sky in Stokes /, Q, and U parameters, convolves this simulated sky with the beam 
pattern for each differencing assembly, then adds uncorrelated instrument noise to each 
pixel in each map. We then co-add the simulated skies in each frequency band and compute 
C I( ^(9) using the same software for both the WMAP data and the si mulations. AH analysi s 



uses only pixels outside the WMAP KpO foreground emission mask ((Bennett et all 12003d ) , 
approximately 76% of the full sky. 

Figure shows C I( ^(8) derived by co-adding the individual correlation functions for 
the frequencies 41, 61, and 94 GHz (Q, V, and W bands) least likely to be affected by Galactic 
foregrounds. The grey band shows the 68% confidence interval for the null simulations. It 
is clear that WMAP detects a temperature-polarization signal at high statistical confidence, 
and that signals exist on both large and small angular scales. We define a goodness-of-fit 
statistic 

X 2 — X! t^MAP ~~ (C S / im)]a M^ 1 [C^j A p — (C,fi m )]&, (7) 
ab 

where C^ AP is the co-added correlation function from WMAP data, (C^) is the mean from 
the Monte Carlo simulations, and M is the covariance matrix between angular bins a and b 
derived from the simulations. We find \ 2 = 207 for 78 degrees of freedom when comparing 
WMAP to the null model: WMAP detects temperature-polarization correlations significant 
at more than 10 standard deviations. 
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Figure 1: Temperature-polarization correlation function for WMAP co-added QVW data. 
The gray band shows the 68% confidence interval for similar co-added data taken from Monte 
Carlo simulations without polarization. The inset shows data for 9 < 10°. The data are 
inconsistent with no temperature-polarization cross-correlations at more than 10 standard 
deviations. Note that the data are not independent between angular bins. 
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2.1 Systematic Error Analysis 

Having detected a significant signal in the data, we must determine whether this signal 
has a cosmological origin or results from systematic errors or foreground sources. We test 
the convergence of the mapping algorithm using end-to-end simulations, comparing maps 
derived from simulated time-ordered data to the input maps used to generate the simulated 
time series. The simulations include all major instrumental effects, including beam ellipticity, 
radiometer performance, and instrument noise (including 1/ f component), and are p rocessed 
using the same map-making software as the WMAP data (jHinshaw et al. . 2003a ). The Q 



and U maps converge rapidly, within the 30 iterations required to derive the calibration 
solution. Correlations in the time-ordered data introduce an anti-correlation in the U map 
at angles corresponding to the beam separation, with amplitude 0.5% of the noise in the map. 
This effect is independent for each radiometer and does not affect temperature-polarization 
cross-correlations. Similarly, residual 1/ / noise in the time series can create faint striping in 
the maps, but does not affect cross-correlations. 

The largest potential systematic error in the temperature-polarization cross-correlation 
results from bandpass mismatches in the amplification/detection chains. We calibrate the 
WMAP data in thermodynamic temperature using the Doppler dip ole from t he satellite's or- 
bit about the Sun as a beam-filling calibration source (jHinshaw et all 2003a ). Astrophysical 



sources with a spectrum other than a 2.7 K blackbody are thus slightly mis-calibrated. The 
amplitude is dependent on the product of the source spectrum with the unique bandpass of 
each radiometer. If the bandpasses in each radiometer were identical, the effect would cancel 
for any frequency spectrum, but differences in the bandpasses between the two radiometers 
in each DA generate a non-zero residual in the difference signal used to generate polarization 
maps (Eq. 0]). This signal is spatially correlated with the unpolarized foreground intensity 
but is independent of the orientation of the radiometers on the sky (polarization angle 7). 
In the limit of uniform sampling of 7 this term drops out of the sky map solution. However, 
the WMAP scan pattern does not view each pixel in all orientations; unpolarized emission 
with a non-CMB spectrum can thus be aliased into polarization if the bandpasses of the two 
radiometers in each DA are not identical. This is a significant problem only at 23 GHz (K 
band), where the foregrounds are brightest and the bandpass mismatch is largest. 

We quantify the effect of bandpass mismatch using end-to-end simulations. For each 
time-ordered sample, we compute the signal in each radiometer using an unpolarized fore - 
ground model and the measured pass bands in each output channel f Jarosik et al. . 2003). 



We then generate maps from the simulated data using the WMAP one-year sky coverage and 
compute C I( ~t{6) using the output /, Q, and U maps from the simulation. We treat this as 
an angular template and compute the least-squares fit of the WMAP data to this bandpass 
template to determine the amplitude of the effect in the observed correlation functions. We 
correct the WMAP correlation functions C I( ^(6) and C IU (6) at K and Ka bands by sub- 
tracting the best-fit template amplitudes. The fitted signal has peak amplitude of 8 fiK 2 at 
23 GHz and 5 fiK 2 at 33 GHz. No other channel has a statistically significant detection of 
this effect. 

Sidelobe pickup of polarized emission from the Galactic plane can also produce spu- 
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rious polarization at high latitudes in the Q and U maps. We esti mate this effec t usin g 
the measured far-sidelobe response for each beam in each polarization ( Barnes et all 12003). 
Sidelobe pickup of polarized structure in the Galactic plane is less than 1 /zK 2 in C I( ^(9) 
at 23 GHz and below 0.1 fj,K 2 in all other bands. We correct the polarization maps for the 
estimated sidelobe signal and propagate the associated systematic uncertainty throughout 
our analysis. Note that all of these systematic errors depend on the Galactic foregrounds, 
and have different frequency dependence than CMB polarization. 

Other instrumental effects are negligible. We measure polarization by differencing 
the outputs of the two radiometers in each differencing assembly (Eq. EJ). Calibration 
errors (as opposed to the bandpass effect discussed above) can alias temperature anisotropy 
into a spurious polarization signal. We have simulated the uncertainty in the calibration 
solution using both rea listic gain drifts and drifts ten times larger than observed in flight 
( Hinshaw et all l2003a| ). Gain drifts (either intrinsic or thermally-induced) contribute less 
than 1 fxK 2 to C IQ (9) in the worst band. 

Null tests provide an additional check for systematic errors. Thomson scattering of 
scalar temperature anisotropy produces a curl-free polarization pattern. A non-zero cosmo- 
logical signal is thus expected only for the IQ (TE) correlation, whereas systematic errors 
or foreground sources can affect both the IQ and IU (TB) correlations. A \ 2 analysis shows 
C IU {9) to be consistent with instrument noise. We further limit possible systematic effects 
by correlating the Stokes I sum map from the Q- or V-band (as noted above) with the 
polarization difference maps (Ql - Q2)/2, (VI - V2)/2, (Wl - W2)/2, and (W3 - W4)/2. 
The temperature (Stokes I) map in all cases is a sum map; the test is thus primarily sensi- 
tive to systematic errors in the polarization data. The difference maps are consistent with 
instrument noise. 



2.2 Foregrounds 

Galactic emission is not a strong contaminant for CMB temperature anisotropy, but could be 
significant in polarization. WMAP measurements of unpolarized foreground emission show 
synchrotron, free-free, and thermal dust emission all sharing significant spatial structure. Of 
these components, only synchrotron emission is expected to generate significant polarization; 
other sources such a s spinn ing dust are limited to less than 5% of the total intensity at 33 



GHz fjBennett et alll2003eft . 

Foreground polarization above 40 GHz is faint: fitting the correlation functions at 
41, 61, and 94 GHz (Q, V, and W bands) to a single power-law C IQ (6, v) = Cq Q (6) {v/v Q f 
yields spectral index (3 = —0.4 ± 0.4, consistent with a CMB signal ((3 = 0) and inconsistent 
with the spectral indices expected for synchrotron ((3 ~ —3), spinning dust (/3 ~ —2), or 
thermal dust (f3 ~ 2). The measured signal can not be produced solely by a single foreground 
emission component. 



A two-component fit 



C^(M=QSb(0) + CgiW(-) (8) 
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tests for the superposition of a CMB component with a single foreground component. Fig- 
ure El shows the resulting decomposition into CMB and foreground components. We obtain 
a marginal detection of foreground component with best-fit spectral index f3 = —3.7 ± 0.8 
consistent with synchrotron emission. We test for consistency or possible residual systematic 
errors by repeating the fit using different temperature maps and different combinations of 
WMAP polarization channels. The fitted CMB component (left panels of Fig. |2J) is robust 
against all combinations of frequency channels and fitting techniques. Note the agreement 
in Fig. |21 between nearly independent data sets: the co-added QVW data (uncorrected for 
foreground emission) and the KKaQ data (corrected for foreground emission). We obtain ad- 
ditional confirmation by replacing the V-band temperature map in the cross-correlation (Eq. 
IHJ) with th e "internal linear com bination" temperature map designed to suppress foreground 



emission ([Bennett et all 12003d ). The fitted CMB component does not change. We test 



for systematic errors by replacing the temperature map with the COBE-DMR map of the 
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Figure 2: Fitted CMB (left) and foreground (right) components from a multi-frequency 
decomposition of the measured two-point correlation functions. Top panels show the IQ 
(TE) correlation, while bottom panels show IU (TB). The CMB component is shown in 
units of thermodynamic temperature, while the foreground is shown in antenna tempera- 
ture evaluated at 41 GHz. Different colors show the effect of using different temperature 
maps in the cross-correlation, or including different polarization frequency channels in the 
CMB- foreground decomposition (see text). The fitted CMB component is stable as different 
frequency channels and data sets are analyzed. Foreground emission is faint compared to 
the cosmic signal. 
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CMB temperature ( Bennett et al. . 19961 ). excluding any instrumental correlation between 



the temperature and polarization data. Again, the results are unchanged. 

We further constrain foreground contributions by computing the cross- correlation 
between the WMAP polarization data and temperature maps dominated by foregrounds. We 
replace the temperature map in Eq. IHlwith either the WMAP maximum-entropy foreground 



model ((Bennett et all 12003d ) or a "residual" foreground map created by subtracting the 



internal linear combination CMB map from the individual WMAP temperature maps. We 
then correlate the foreground temperature map against the WMAP polarization data in 
each frequency band, and fit the resulting correlation functions to CMB and foreground 
components (Eq. IBJ). The two foreground maps provide nearly identical results. The fitted 
CMB component has nearly zero amplitude, consistent with the instrument noise. The 
fitted foreground has amplitude 0.5 ± 0.1 /iK 2 at Uq = 41 GHz, with best-fit index (3 = —3.4 
consistent with synchrotron emission. 



3 POLARIZATION CROSS-POWER SPECTRA 



In a second analysis method, we compute the angular power spectru m of the temperatur e- 
polarization correlations using a quadratic estimator (cf Appendix A in Kogut et al. ( 2003h ). 
We compute cf E and cJ B individually for the each WMAP frequency band, using uniform 
weight for the temperature map and noise weight for the polarization maps. We then combine 
the angular power spectra, using noise- weighted QVW data for / > 21 where foregrounds 
are insignificant, and a fit to CMB plus foregrounds using all 5 frequency bands for Z 5 21. 
Since foreground contamination is weak, we gain additional sensitivity in this analysis by 
using the Kp2 sky cut retaining 85% of the sky. 

We estimate the uncertainty in each I bin using the covariance matrix M for the 
polarization cross-powe r spectrum. Based on our analysis of the c[ T covariance matrix 



(iHinshaw et all l2003bh . ;he cf covariance matrix has the form along the diagonal of 



M, 



< c; 



TETE 
I c l 



> 



< C 



TE 



> 



(cJ T + n TT / Wl )(cf E + Uee/wi) + (cJ E f 



(21 



^■)fskyfsky 



(9) 
(10) 



where utt and uee are t he TT and EE nois e bias terms, wi is the effective window function 



for the combined maps ([Page et all l2003al ). cf T and cf E are the temperature and polar- 



0.85 is the fractional sky coverage for the Kp2 mask, 



and = fsj.y/1.14: fo r noise weighting. We take the c[ T term from the measured tern- 



ization angular power spectra, f s k y 

. :htin g. 

perature power spectra (IHinshaw et all |2003bJ) and the cf E term predicted by the best-fit 
ACDM model (|Soergel et all l2003h fallowing cf E to vary as a function of optical depth in 
the likelihood analysis). Monte Carlo simulations demonstrate that the analytic expression 
accurately describes the diagonal elements of the covariance matrix. We approximate the 
off-diagonal terms using the geometric mean of the covariance matrix terms for uniform and 
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noise weighting ((Hinshaw et all l2003bl ). 3 



r Al- 



The largest off-diagonal contribution, —2.8%, is at Al = 2 from the symmetry of our sky 
cut and noise coverage. The total anticorrelation is Z^A/^o r A« = —0.124. Because of this 
anti-correlation, the error bars for the binned c[ E are slightly smaller than the naive estimate. 

Figure 01 shows the polarization cross-power spectra for the WMAP one-year data. 
The solid line shows the predicted signal for adiabatic CMB pertu rbations, based onl 



on a fit to the measu red temperature angular power spectrum ci (jSpergel et all 1200 



Hinshaw et al 



Two features are apparent. The TE data on degree a ngular scales 



(7 > 20) are in excellent agreement with a priori predictions of adiabatic models (jCoulson et al 
1994 ) . Other than the specification of adiabatic perturbations, there are no free parame- 

of 24.2 for 23 degrees of freedom indicates 



ters 



the solid line is not a fit to c[ E . The x 2 



that the CMB anisotropy is dominated by adiabatic perturbations. On large angular scales 
(/ < 20) the data show excess power compared to adiabatic models, suggesting significant 
reionization. 

The WMAP detection of the acoustic structure in the TE spectrum confirms several 
basic elements of the standard paradigm. The amplitudes of the peak and anti-peak are a 
measure of the thickness of the decoupling surface, while the shape confirms the assumption 

3 Note that Hinshaw ct al. ( 2003b) define off-diagonal elements in terms of the inverse covariance matrix, 
which differs from r&i by a sign. 
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Figure 3: Polarization cross-power spectra cj for the WMAP one-year data. Note that we 
plot (I + l)/27r cf E and not 1(1 + l)/27r cf E . This choice emphasizes the oscillatory nature 
of cJ E . For clarity, the dotted line shows q = 0. The solid line is the predicted signal based 
on the cJ T power spectrum of temperature anisotropy - there are no free parameters. The 
TE correlation on degree angular scales (7 > 20) is in excellent agreement with the signal 
expected from adiabatic CMB perturbations. The excess power at low / indicates significant 
reionization at large angular scales. 
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that the primordial fluctuations are adiabatic. Adiabatic fluctuations predict a tempera- 
ture/polarization signal anticorrelated on lar ge scales, with TE pea ks and anti-peaks lo- 
cated midway between the temperature peaks iHu fc Sugivamal (|l994| ). The existence of TE 
correlations on degree angular scales also provides evidence for super-horizon temperature 
fluctu ations at decoupling, as expected for inflationary models of cosmology (jPeiris et al 
2003h 



4 REIONIZATION 

WMAP detects statistically significant correlations between the CMB temperature and po- 
larization. The signal on degree angular scales (I > 20) agrees with the signal expected in 
adiabatic models based solely on the temperature power spectrum, without any additional 
free parameters. We also detect power on large angular scales (I < 10) well in excess of the 
signal predicted by the temperature power spectrum alone. This signal can not be explained 
by data processing, systematic errors, or foreground polarization, and has a frequency spec- 
trum consistent with a cosmological origin. 

The signal on large angular scales has a natural interpretation as the signature of 
early reionization. Both the temperature and temperature-polarization power spectra ca n 



be related to the power spectrum of the radiation field during scattering ( Za ldarriag a. 1997). 
Thomson scattering damps the temperature anisotropy and regenerates a polarized signal 
on scales comparable to the horizon. The existence of polarization on scales much larger 
than the acoustic horizon at decoupling implies significant scattering at more recent epochs. 

4.1 Reionization in a ACDM Universe 

If we assume that the ACDM model is the best description of the physics of the early universe, 
we can fit the observed temperature-polarization cross-power spectrum to derive the optical 
depth r. We assume a step func tion for the ionization fraction x e and use the CMBFAST 



^ep i 

code ((Seljak fc Zaldarriagal 199fi[ ) to predict the multipole moments as a function of optical 
depth. While this assumption is simplistic, our conclusions on optical depth are not very 
sensitive to details of the reionization history or the background cosmology. 

Figure0]compares the polarization cross-power spectrum c[ E derived from the quadratic 
estimator to ACDM models with and without reionization. The rise in power for Z < 10 is 
clearly inconsistent with no reionization. We quantify this using a maximum-likelihood 
analysis 

e*p(-k 2 ) , 12) 

| M |l/2 • ^ 

Figure El shows the relative likelihood £/Max(£) for the optical depth r assuming a ACDM 
cosmology, with all other parameters fi xed at the values derived from the temperature power 



spectrum alone (jSpergel et all |2£)03). The likelihood for the 5-band data corrected for 
foreground emission peaks at r = 0.17 ± 0.03 (statistical error only): WMAP detects the 
signal from reionization at high statistical confidence. 
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A full error analysis for r must account for systematic errors and foreground un- 
certainties. We propagate these effects by repeating the maximum likelihood analysis using 
different combinations of WMAP frequency bands and different systematic error corrections. 
We correct C I( ^(9) in each frequency band not for the best estimate of the systematic error 
templates, but rather the best estimate plus or minus one standard deviation. We then fit 
the mis-corrected C I( ^(9,u) for a CMB piece plus a foreground piece (Eq. |HJ) and use the 
CMB piece in a maximum-likelihood analysis for r. The change in the best-fit value for r 
as we vary the systematic error corrections propagates the uncertainties in these corrections. 
Systematic errors have a negligible effect on the fitted optical depth; altering the systematic 
error corrections changes the best-fit values of r by less than 0.01. 

The largest non-random uncertainty is the foreground separation. We assess the 
uncertainty in the foreground separation by repeating the entire systematic error analysis 
(using both standard and altered systematic error corrections) with the foreground spectral 
index (3 = — 3.7±0.8 shifted one standard deviation up or down from the best-fit value. Fitted 
values for r derived from the three high-frequency channels (QVW) without foreground 
fitting are in agreement with lower- frequency data once foregrounds are taken into account. 
We obtain nearly identical values for r when fitting either the highest-frequency data set 
QVW or the lowest-frequency set KKaQ. The fitted optical depth is insensitive to the spectral 
index: varying the spectral index from -2.9 to -4.5 changes the fitted values by 0.02 or less. 
We adopt r = 0.17 ± 0.04 as the best estimate for the optical depth to reionization, where 
the error bar reflects a 68% confidence level interval including statistical, systematic, and 
foreground uncertainties. 
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Figure 4: WMAP Polarization cross-power spectra cJ E (filled circles) compared to ACDM 
models with and without reionization. The rise in power for I < 10 is consistent with reion- 
ization optical depth r = 0.17 ± 0.04. The error bars on WMAP data reflect measurement 
errors only; adjacent points are slightly anti-correlated. The grey band shows the 68% con- 
fidence interval from cosmic variance. The value at I — 7 is particularly sensitive to the 
foreground correction. 
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Soergel et all (|2003l) include the TE data in a maximum-likelihood analysis combining 



WMAP data with other astronomical measurements. The resulting value, r = 0.17 ±0.06, is 
consistent with the value derived from the TE data alone. The larger uncertainty reflects the 
effect of simultaneously fitting multiple parameters. The TE analysis propagates foreground 
uncertainties by re-evaluating the likelihood using different foreground spectral index. Since 
foreground affect only the lowest multipoles, the combined analysis propagates foreground 
uncertainty by doubling the statistical uncertainty in cf E for 2< IliA to account for this effect. 



4.2 Model-Independent Estimate 

An alternative approach avoids assuming any cosmological model and uses the measured 
temperature angular correlation function to determine the radiation power spectrum at re- 
combination. This approach assumes that the best estimate of the three dimensional radia- 
tion power spectrum is the measured angular power spectrum rather than a model fit to the 
angular power spectrum. Given the observed temperature power spectrum cj T , we derive the 
predicted polarization cross-power spectrum c[ E , which we then fit to the observed TE spec- 
trum as a function of optical depth r. We obtain r = 0.16 ±0.04, in excellent agreement with 
the value derived assuming a ACDM cosmology. We emphasize that the model-independent 
technique makes no assumptions about the cosmology. The fact that it agrees well w i th the 
best-fit model from the combined temperature and polarization data (jSpergel et all l200ah 
is an additional indication that the observed temperature-polarization correlations on large 
angular scales represent the imprint of physical conditions at reionization. The dependence 
on the underlying cosmology is small. 
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Figure 5: Likelihood function for optical depth r for a ACDM cosmology, using all 5 WMAP 
frequency bands fitted to CMB plus foregrounds with foreground spectral index (3 = —3.7. 
After including systematic and foreground uncertainties the optical depth is consistent with 
a value r = 0.17 with 95% confidence range 0.09 = rl 0.28. 
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4.3 Early Star Formation 



Reionization can also be expressed as a redshift z r assuming an ionization history. We 
consider two simple cases. For instantaneous reionization with ionization fraction x e — 1 
at z < z r , the measured optical depth corresponds to redshift z r = 17 ± 3. This conflicts 
with measurements of the Gunn-Peterson absorp tion trough i n spectra of distant quasars, 
which show neutral hydrogen present at z ~ 6 (|Beckerl . 2001 : Diorgovski et al. . 2001 : Fan . 
2002:[). Reionization clearly did not occur through a single rapid phase transition. However, 
since absorption spectra are sensitive to even small amounts of neutral hydrogen, models 
with partial ionization have enough neutral column density to produce the Gunn- 

Peterson trough while still providing free electrons to scatter CMB photons and produce 
large-scale pola rization. D irect Gunn-Peterson observations only imply a neutral hydrogen 



fraction > 1% (Fa 



ify the simplest model to add a second 
transition: a jump from x e = to x e = 0.5 at redshift z r , followed by a second transition 
from x e = 0.5 to x e — 1 at redshift z — 7. Fitting this model to the measured optical depth 
yields z r « 20. In reality, reionization is more complicated than simple step transitions. 
Allowing for model uncertainty, the measured optical depth is consistent with reionization 
at redshift 11 < z r < 30, corresponding to times 100 < t r < 400 Myr after the Big Bang 
(95% confidence). 

Extrapolations of the observe d ionizing flux to highe r redshift lead to predicted CMB 
optical depth between 0.04 — 0.08 ( Miralda-Escudel . l2002[ ). lower than our best fit values. 
The measured optical depth thus implies additional sources of ionizing flux at high redshift. 
An early generation of v ery massive (Pop III) stars could provide the required additional 
heating. iTegmarkl (jl997h estimate that 10~ 3 of all baryons should be in collapsed objects 
by z = 30. If these baryons form massive stars, they would reionize the universe. How- 
ever, photons below the hydrogen ionization threshold will destroy molecular hydrogen (the 
principal vehicle for cooling in early stars), driving the effective mass t hreshold for star for- 
mation to ~ 10 8 sol ar masses and impedin g subsequent star formation (lHaiman et al.Lll997 : 



Gnedin Sz Ostriker . 1997 : Tegmarkl Il997t ). X-ray heating and ionization ( Venkatesan et al._ 
200lt lOhl 2001 ) may provide a loophole to this argument by enhancing the formation of H 2 



molec ules (lHaiman et al. . 2000i ). 



CeJ (|2n03) provides a physically-motivated model of "double reionization" that re- 
sembles the two-step model above. A first generation of massive Pop III stars initially 
ionizes the intergalactic medium. The increased metallicity of the intergalactic medium then 
produces a transition to smaller Pop II stars, after which the reduced ionizing flux allows 
regeneration of a neutral hydrogen fraction. The ionization fraction remains at x e ~ 0.6 
until the global star formation rate surpasses the recombination rate at z — 6, restoring 
x e = 1. The predicted value r = 0.10 ± 0.03 should be increased somewhat to reflect the 



higher WMAP values for the baryon density Qf, and normalization o~s (Spcrgel et al., 2003 



The contribution from ionized helium will also serve to increase r (Venkatesan et al., 200 



)03; 



Wvithe fc Loebl |2003| ). The WMAP determination of the optical depth indicates that ion- 
ization history must be more complicated than a simple instantaneous step function. While 
physically plausible models can reproduce the observed optical depth, reionization remains 
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a complex process and can not be fully characterized by a single number. A more complete 
determi nation of the ionization history requires evalu ation of the detailed TE and EE power 
spectra (jKaplinghat et all Eq3 IHu fc Holder! HT" 



The WMAP mission is made possible by the support of the Office of Space Sciences 
at NASA Headquarters and by the hard and capable work of scores of scientists, engineers, 
technicians, machinists, data analysts, budget analysts, managers, administrative staff, and 
reviewers. 
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